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PACS 97. lO.Gz - Accretion and accretion disks stellar 
PACS 95 . 30 . Qd - Magnetohydrodynamics in astrophysics 

Abstract - Accretion disks around very compact objects such as very massive Black hole can 
grow according to thick toroidal models. We face the problem of defining how does change 
the thickness of a toroidal accretion disk spinning around a Schwarzschild Black hole under the 
influence of a toroidal magnetic field and by varying the fluid angular momentum. We consider 
both an hydrodynamic and a magnetohydrodynamic disk based on the Polish doughnut thick 
model. We show that the torus thickness remains basically unaffected but tends to increase or 
decrease slightly depending on the balance of the magnetic, gravitational and centrifugal effects 
which the disk is subjected to. 
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Introduction. — Configurations of extended matter 
rotating around gravitational sources are undoubtedly an 
environment with a rich phenomenology covering very dif- 
ferent aspects of Astrophysics in different energetic sec- 
tors, from for example the proto-planetary accretion disks 
to, most likely, violent dynamical effects such as Gamma 
Ray Bursts (GRB), see e.g. [1,2]. In addition, since many 
of these disks are made of plasma and they share their 
toroidal topology with other plasma devices in the terres- 
trial physic laboratory , for example Tokamak machines, 
the questions about the plasma stability and the confine- 
ment under magnetic field, involve in general a much wider 
interest of the astrophysic phenomenology. There are dif- 
ferent ways to classify accretion disk models, for example 
according to optical depth, accretion rates and geometric 
structure: in general accretion disks have toroidal shape, 
but we can roughly distinguish between thick and thin 
disk models, in this work we consider a thick toroidal 
disk. The dynamics of any accretion disk is determined 
by several factors, say centrifugal, dissipative, magnetic 
effects. In a sense we could say that the thick accretion 
disks are strongly characterized by the gravitational ones 
in the following respects: thick models are supposed to 
require attractors that are very compact sources of strong 
gravitational fields, for example Black hole and they char- 
acterize the physics of most energetic astrophysical objects 
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as Active Galactic Nuclei (AGN) or GRB. Since they char- 
acterize regions very close to the source, it is in many cases 
necessary a full general relativistic treatment of the model. 
Furthermore in many thick models gravitational force con- 
stitutes the basic ingredient of the accretion mechanism 
independently of any dissipative effects that are strate- 
gically important for the accretion process in the thin 
models [3,4]. The Polish doughnut (PD) is an opaque 
(large optical depth) and Super-Eddington (high matter 
accretion rates) disk, characterized by an ad hoc distri- 
butions of constant angular momentum. This model and 
its derivations are widely studied in the literature with 
both numerical and analytical methods, we refer for an 
extensive bibliography to [5,6]. The torus shape is defined 
by the constant Boyer potential W: closed equipotential 
surfaces define stationary equilibrium configurations, the 
fluid can fill any closed surface. The open equipotential 
surfaces define dynamical situations, for example the for- 
mation of matter jets [7]. The critical, self-crossing and 
closed equipotential surfaces (cusp) locates the accretion 
onto the Black hole due to Paczynski mechanism: a little 
overcoming of the critical equipotential surface, violation 
of the hydrostatic equilibrium when the disk surface ex- 
ceeds the critical equipotential surface Wcusp- The rela- 
tivistic Roche lobe overflow at the cusp of the equipoten- 
tial surfaces is also the stabilizing mechanism against the 
thermal and viscous instabilities locally, and against the so 
called Papaloizou and Pringle instability globally [8] . In [9] 
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and [10], the fully relativistic theory of stationary axisym- 
metric PD torus in Kerr metric has been generalized by 
including a strong toroidal magnetic field, leading to the 
analytic solutions for barotropic tori with constant angu- 
lar momentum. In a recent work the PD model has been 
revisited providing a comprehensive analytical description 
of the PD structure, considering a perfect fluid circularly 
orbiting around a Schwarzschild background geometry, by 
the effective potential approach for the exact gravitational 
and centrifugal effects [6]. The work has taken massively 
advantage of the structural analogy of the potential for 
the fluid with the corresponding function for the free par- 
ticle dynamics which is not subjected to kinetic pressure 
in the same orbital configuration around the source. This 
analysis leads to a detailed, analytical description of the 
accretion disk, its toroidal surface, the thickness, the dis- 
tance from the source with the effective potential and the 
fluid angular momentum. A series of results concerning 
the PD shape was found, in particular it was proved that 
the torus thickness {h) increases with the "energy" pa- 
rameter (K), defined as the value of the effective poten- 
tial for that momentum, and decreases with the angular 
momentum: the torus becomes thinner for high angular 
momenta, and thicker for high energies and the maximum 
diameter (A) increases with K, but decreases with the fiuid 
angular momentum, moreover the location of maximum 
thickness of the torus moves towards the external regions 
with increasing angular momentum and energy, until it 
reaches a maximum an then decreases. 

In this paper we generalize this approach to the PD disk 
to the magnetohydrodynamic models developed in [9,10]. 
We address the specific question of torus squeezing on the 
equatorial plane exploring the disk thickness changing the 
physical characteristics of the torus as the angular mo- 
mentum and the effective potential which characterize the 
toroidal surface, and the magnetic field for the magnetohy- 
drodynamic case, then we will find the extreme squeezing. 

Fluid configuration. — We consider the motion in 
the Schwarzschild background geometry written as 

ds^ = -e'^^^dt' + e-^^^\r)dr'^ + [dO'^ + sin' Od^^) , 

(1) 

in the spherical polar coordinate (t, r, ^, (j)) where e^^^^ = 
(l-2M/r). We define A = /7^ S = t/S ^ = /7^, 
6 = [/^, for the fiuid four- velocity vector field [/, and 
we introduce the set of variables {E^ Vsc-, T} by the fol- 
lowing relations^: 



E 



is the effective potential [6] . We consider the case of a fiuid 
circular configuration, defined by the constraints A = 
(i.e. Vsc = E)^ restricted to a fixed plane sin6> = a ^ 0. 
No motion is assumed in the angular direction, which 
means = 0. For the symmetries of the problem, we 
always assume 9tQ = and S^Q = 0, being Q a generic 
tensor of the spacetime. The fiuid is characterized by the 
constant angular momentum: 
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We consider an infinitely conductive plasma: FabU^ = 
0, where Fab is the Faraday tensor. Using the equation 
U^Ba = 0, where is the magnetic field, we obtain the 
relation B^ = <I>5^/I]. We note that, if we set 5^ = 
from the Maxwell equations, we infer B^ cot ^ = (with 
£ =constant), that is satisfied for 5^ = or ^ = 7r/2. 
The former implies that the assumptions dc/^B^ = and 
B^ = together lead to 5^ = 0, in this work we consider 
90^^ = and B^ = B^ = 0. As noted in [9] the presence 
of a magnetic field with a predominant toroidal component 
is reduced to the disk differential rotation, viewed as a 
generating mechanism of the magnetic field, for further 
discussion we refer to [9-11] and [12-15]. 

The Euler equation for this system has been exactly 
integrated for the background spacetime of Schwarzschild 
and Kerr Black holes in [9-11] assuming the magnetic field 
is 
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where ps = M (gtc/^gtc/) - gttgc/>c/>)'^~^^'^ is the magnetic 
pressure, uj is the fiuid enthalpy, q and A4 are constant, 
we assume moreover a barotropic equation of state. We 
can write the Euler equation as 



drP 



dep 



P^P p^p 
where using eq. (4), it follows 
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= InVsc, Wy^^ =gr{r,e)+gr{0), (8) 
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where gei^^) and gr{0) are functions to be fixed by the 
(3) integration. For ^ 7^ 1 it is Grir^O) = Geir^O) = G{r,e), 
where 



^We adopt the geometrical units c = G = M = 1 extended to the 
electromagnetic quantities by setting eo = 1, where eo is the vacuum 
permittivity and the (— , +, +, +) signature. 
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Thus, the general integral is, in the magnetic case and for 

/ J^ = _(Vr(/)+Vr(e-)), (11) 
J P + P 



Wr = g{r,e) + \n{V,e)+9r{0), 
We = g{r,e) + \n{V,e)+9e{r). 



(12) 



The Boyer surfaces. — The procedure described in 
the present article borrows from the Boyer theory on the 
equipressure surfaces for a PD torus [6,7], we consider the 
equation for the Boyer potential: W = Q{r^O) -\- In^l^c) = 
c = constant, where degr{0) = drge{r) and setting = 
Qr = 0, we can write 



Ke- 



(13) 



where K = e^, the toroidal surfaces are obtained from 
the equipotential surfaces [6,7]. eq. (13) can be solved 
for the radius r assuming q ^ 0. We notice that, in the 
limit case ^ = 0, the magnetic field 5, does not depend 
on the fluid enthalpy, furthermore eq. (13) in this case is 
Vsc = K, this means that the magnetic fleld \ q^Q does 
not affect in anyway the Boyer potential and therefore the 
Boyer surfaces. 

Consider now the case of small eq. (13) is 
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assuming the enthalpy a; to be a constant, we introduce 
the following x = rcos6>, y = rsin6>, and S = Mq/uo. 
In the following we consider the approximation <S ~ 0, we 
note that the ratio M./uo gives the comparison between the 
magnetic contribution to the fluid dynamics, through A^, 
and the hydrodynamic contribution through its speciflc 
enthalpy uo. Solving eq. (14), written in coordinates (x, ?/), 
for X we obtain the Boyer surfaces: 
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The coordinates of the maximum point of the surface sat- 
isfy the condition dx/dy = and d^x/d^y < 0, this equa- 
tion can be solved in order to obtain the exact form of 
{xmax^ymax)^ the torus thickucss is here deflned as the 
maximum height of the surface i.e. h = 2xmax with coor- 
dinate ymax' We then deflne the ratio 



where Rs{S = 0) is the squeezing in the PD case {S = 0), 
one can see this term as the hydrodynamic contribution 
to the disk squeezing. Where the approximation 5 <C 1 
holds we can study the linear dependence of Rs from S. 
The case of increasing (decreasing ) squeezing with S is 
deflned by (^ 1(^=0) ) > (^ 1(^=0) < 0). The Boyer 
surfaces for 5 = and 5 7^ are shown in flgs. 1. 

Remark: The magnetic pressure is regarded as a per- 
turbation of the hydrodynamic component, it is assumed 
that the Boyer theory remains valid and applicable in this 
approximation (discussions on similar assumptions can be 
found in [5, 9]): the torus shape is determined by the 
equipotential surfaces which now are regulated in eq. (12) 
by an effective potential deformed, respect to the PD case 
(5 = 0), by the magnetic fleld. The approximation 5 ~ 
follows from this assumption. This analysis should be seen 
as a flrst attempt to compare the magnetic contribution to 
the torus dynamics respect to the kinetic pressure in terms 
of the torus squeezing. However, we have performed sev- 
eral numerical analysis by setting a solution of the family 
of magnetic flelds assigning a flxed value for the parameter 
q. The range of this parameter is clearly divided into two 
regions, say < q < 1 and g > 1, with a subrange with 
extreme g = 2. The case q = 2 is indeed interesting be- 
cause the magnetic fleld loops wrap around with toroidal 
topology along the torus surface. 

For q < 1 the trend is clear: closed surfaces are roughly 
those studied in the case ^ 0, for > 1 we get closed 
surfaces in the limit 5 ~ 0, which is in agreement with 
expansion around 5 = 0, that is, when the contribution 
of the magnetic pressure to the torus dynamics is prop- 
erly regarded as a perturbation with respect to the hydro- 
dynamics solution, this is realized when an appropriate 
condition (5 <C 1) on the fleld parameters is satisfled. 
First we introduced the parameter S = MquJ^~^/{q — 1) 
in eq. (13), then we expand eq. (13) around 5 = 0, and 
collect the terms up to flrst order. It should be noted that 
this equation, is equivalent to the corresponding obtained 
by an expansion in q = 0, once ignored the terms KSq, and 
changing S —S (by convention we are taking S > 0). 
The system analysed in eq. (15) is thus equivalent to the 



case of a torus with S = 



Rs=y 



(16) 



as the squeezing function for the torus, where X = ys — y2 
is the maximum diameter of the torus surface section, see 
flg. 1, and (^1,^2,^3), are solutions of the equations x'^ = 
0. 

Consider i^^, for 5 <C 1: 



^Oand {qKS) = 0[S] 
(i.e. higher order than the 1th degree). 

Squeezing of the Polish doughnuts. — In this sec- 
tion we consider the squeezing of the PD model {B^ = 0). 
The PD torus has been extensively studied in the litera- 
ture with both analytical and numerical methods, we refer 
in particular to [6] with which we share much of the nota- 
tion and conventions used in this work. Setting 5 = in 
eq. (15) we flnd 



2{KH^^y^) 

K2 (£2_^2)^^2 
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Rs ^ Rs 1(5=0) + 1(5=0) 



(17) 



with 2^/2/3 < K < 1 and angular momentum (^;. )^ < 
i'^ < (^^)^ where > 27/2, we refer to [6] for the exact 
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K=0.98 /^=18 K=0.981 f=\6 
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Fig. 1: (Colour on-line)Left Case — 0: the closed Boyer surface at K = 0.981 and = 18: X = ys — y2 is the surface 
maximum diameter, ^ = ^2 — 2 is the distance from the source, 6 = y2 — yi is the distance from the inner surface, 6 = yi — 2 is 
the distance of the inner surface from the horizon, h = 2xm is the maximum thickness of the torus (torus height). Right Case 
/ 0: torus surface at K = 0.981 and i'^ — 16, in the plane (r sin^, r cos^) for different values of S. Inset figures are enlarged 
view of the main plot. 



expression of and details, fig. 1-right shows a section of 
the PD torus for K = 0.981 and = 18. We summarize 
the results concerning the squeezing function Rg = h/X'm. 
the following. 

The squeezing function Rg increases monotonically with 
K and, at fixed K decreases with figs. 2. This evolution 
can also be seen by studying of the function Ag = h — X, 
shown in figs. 3,4 as function of £ and K. However we 
note that the lower Rg = h/X is and the thinner is the 
torus, conversely higher is Rg and thicker is the torus. 
This therefore means that the toroidal surface is squeezed 
on the equatorial plane with decreasing K and increasing 
i'^. These results are indeed consistent with the analysis 
in [6] where it was shown that the torus thickness /i, in- 
creases with the K and decreases with the angular momen- 
tum, that is the torus becomes thinner for high angular 
momenta, and thicker for high energies. The maximum 
diameter A increases with the K, but decreases with the 
fluid angular momentum, moreover the location of max- 
imum thickness of the torus moves towards the external 
regions with increasing angular momentum and K, until 
it reaches a maximum an then decreases. 

We note a region of i and K values, in which this trend 
is reversed, that is Rg = h/X < 1 figs. 2- (left-center), and 
Ag = h — X < figs. 3,4-left, i.e. the squeezing is such 
that the disk is longer that thick, although the minimum 
value Rg ~ 0.95 is still quite large because it still makes 
sense that the model of PD be a thick disk. One sees from 
figs. 4, in this region, a minimum for Ag as a function of i. 
Following the curve for increasing values of K, the function 
Ag decreases or decreases until it reaches a minimum and 
then increases. From fig. 2- center one also sees that the 
lower the angular momentum is and more configurations 
with Rg < 1 are. Finally, some curves overlap in this 
region for one or more values of K, this means that torus 
configurations would exist with different values of but 
equal K and Rg. 

Fig. 2-right shows the solutions of equation i^^ = 1, on 
the curve ^^(K) . Essentially this curve confirms the pre- 
vious results showing that the torus where the {h) equals 



the length on the equatorial plane A increases in angular 
momentum with K, first slowly and then more and more 
rapid as K approaches the value limit K = 1. A fixed Rg, 
the centrifugal effects must compensate for an increase of 
K which tends to bring the disk to a greater Rg, and to 
sufficiently small values of K G (0.97367, 0.9746), the curve 
bends back on itself, see Figs. 2. 

Squeezing of the magnetized torus. — In this sec- 
tion we summarize the analysis concerning the squeezing 
of the magnetized torus {S 7^ 0). 

Fig. 4 shows several surfaces Rg{£'^,K) for different val- 
ues of an enlarged view of this plot, distinguishes the 
different <S-surfaces. In general the function Rg{i'^,K), 
as for the case 5 = 0, increases with K but decreases 
with that is the torus is thicker as the K parameter in- 
creases and becomes thinner as fluid angular momentum 
increases, fig. 4-center, furthermore Rg increases with 5, 
fig. 4,5-left. 

The derivative {dRg/dS){S = 0) is overwhelmingly pos- 
itive: for the magnetized torus, Rg grows with 5, hence the 
torus is thicker as S increases and thinner as S decreases, 
figs. 5,6-left. However, as evidenced by the Figs, 5,6 there 
is a small region 15.3 < i'^ ^ 16.65, in which the deriva- 
tive {dRg/dS){S = 0) is negative and Rg is a decreasing 
function of S: from this we conclude that Rg decreases 
compared to the case <S = i.e. the torus becomes thin- 
ner with increasing of S and viceversa becomes thicker to 
the decrease of 5, figs. 5. These results are confirmed in 
fig. 6, where the squeezing has been plotted as a function 
of S and for different values of K, however it should 
be noted that although the figures show the behavior of 
the ratio Rg across the range < <S < 1, the approxima- 
tions assumed in this work are valuable for small S. Fig. 6 
shows that Rg is generally increasing with K and decreas- 
ing with and also shows the existence of a region where 
this trend is reversed. 

Conclusions. — We have studied the Polish dough- 
nut (PD) model of thick disk for a perfect fluid, circu- 
larly orbiting around a Black hole of Schwarzschild, in 
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Fig. 2: (Colour on-line) Case = 0. Left: The squeezing Rs as function of K and the angular momentum P. The black plane is Rs = 1, 
the curves (^J)^ are plotted. Center: Rg curves for fixes values of the i'^ as function of K, a zoom of the curves is in the inset plot: Rs 
increases monotonically with K and, at fixed K, decreases with Right: curve Rg = 1 (black curve), orange dashed line is (^^)^ black 
dashed line is (^^)^, inside plot is a zoom of the i^s = 1 in the region K G [0.9735, 0.975] 
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Fig. 3: (Colour on-line) Case = 0. Left: difference As = /i — A as function of K and Center plot shows As and the planes 
As = (0, ±2). Left: curves As =cons as function of i"^ and K. 




Fig. 4: (Colour on-line) Left: Case B^ = 0, plot of As as function of K, for selected values of i. Case B^ / 0. Center: the squeezing Rg 
as function of K and the angular momentum with i"^, for selected values of the *S parameter. The black plane is Rs = 1. A zoom is shown 
in right plot: white surface *S = 0, shaded light gray surface *S = 1/8, gray surface S = 1/4, shaded white surface S = 1/2, shaded gray 
surface S = 3/4. 
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Fig. 5: Case B^ / 0. Left: the squeezing surface Rs =cons for selected values of S as function of i'^ and K. The inset plot shows the 
curves Rs = 1 in the region K G (.9765, .98) P G (20, 24). Center: the function R'^(0) = dRs/dS in <S = as function of the P and the 
parameter K. The function is potted in the range R'd^) G (—0.02, +0.02), black plane is R'ciS^) = 0. Right plot is a zoom in the region 
i?^(0) G (-0.02, 0)0, K G (0.954, 1), P G (15.3, 16.65). 
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Fig. 6: Case / 0. Left: curves Rc(0) ^constant in the plane (-^^, K), inside plot is a zoom in the region i'^ G (15.3, 16) where Rc(0) < 0. 
Plot of the squeezing Rg as function of *S for selected values of i'^ and K = 0.97 (center plot), K = 0.99 (right plot). Plots show the range 
*S G [0, 1], however the approximation *S << 1 holds for S < 0.3. 



an hydro dynamic and magnet ohydro dynamic context, us- 
ing the effective potential approach examined in detail for 
the hydrodynamical model in [6]. For the magnetohydro- 
dynamic model we used the solutions of the barotropic 
torus with a toroidal magnetic field discussed in [9, 10]. 
We have considered the magnetic contribution to the to- 
tal pressure of the fluid as a perturbation of the effective 
potential for the exact full general relativist ic part. Once 
defined the torus squeezing function Rg = h/X diS the ratio 
maximum torus height-to-maximum diameter of the torus 
section, we studied Rg varying of the angular momentum 
the effective potential K and the parameter 5, the last 
evaluates the magnetic contribution to the fluid dynam- 
ics respect to the hydrodynamic contribution through its 
specific enthalpy. The lower is Rs and the thinner is the 
torus, conversely the higher is Rs and thicker is the torus. 
It is shown that given the presence of a magnetic pressure 
as a perturbation of the hydrodynamic component, there 
are not quantitatively significant effects on the thick disk 
model; however the analysis showed that the squeezing 
function has in general a monotonic trend with (^, K,tS). 
For the hydrodynamic case the squeezing function Rs in- 
creases monotonically with K, and at fixed K decreases 
with this therefore means that the toroidal surface is 
squeezed on the equatorial plane with decreasing "energy" 
K and increasing these results are indeed in agree- 
ment with analysis in [6]. However we observed a region 
of £ and K values, in which this trend is reversed, that is 
Rs = h/X < 1 reaching a minimum values of Rs ~ 0.95, 
for increasing values of K, the squeezing decreases or de- 
creases until it reaches a minimum and then increases. In 
the magnet ohydrodynamic case, we evaluated the infiu- 
ence of the magnetic field through the parameter S. In 
general, as for the case 5 = 0, the torus is thicker as the 
K parameter increases, and becomes thinner as fiuid an- 
gular momentum increases, furthermore Rs increases with 
S. However there is a small region of ^, in which Rs is a 
decreasing function of S: that is Rs decreases compared to 
the case 5 = i.e. the torus becomes thinner with increas- 
ing of S and viceversa becomes thicker with decreasing S. 

In conclusion, this analysis shows that the torus thick- 
ness is not modified in a quantitatively significant way by 



the presence of a toroidal magnetic field although it is 
much affected by the variation of the angular momentum. 
However we consider this analysis as a first attempt to 
compare the magnetic contribution to the torus dynamics 
with respect to the kinetic pressure in terms of the torus 
squeezing and, more in general, the plasma confinement 
in the disk toroidal surface. 
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